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We show that, within the Quasiparticle Random Phase Approximation (QRPA) and the renormal- 
ized QRPA (RQRPA) based on the Bonn CD nucleon-nucleon interaction, the competition between 
the pairing and the neutron-proton particle-particle and particle-hole interactions causes contribu- 
tions to the neutrinoless double-beta decay matrix element to nearly vanish at internucleon distances 
of more than 2 or 3 fermis. As a result, the matrix element is more sensitive to short-range/high- 
momentum physics than one naively expects. We analyze various ways of treating that physics 
and quantify the uncertainty it produces in the matrix elements, with three different treatments of 
short-range correlations. 

PACS numbers: 23.10.-s; 21.60.-n; 23.40.Bw; 23. 40. He 

Keywords: Neutrino mass; Neutrinoless double beta decay; Nuclear matrix element; Quasiparticle random 
phase approximation 



I. INTRODUCTION 

Neutrino oscillations are firmly established (see, e.g. 
P) S S 0) Hi) and demonstrate that neutrinos have 
masses many orders of magnitude smaller than those 
of charged leptons. But since the masses are nonzero, 
neutrinoless double-beta (Oz^/3/3) decay experiments will 
likely tell us sooner or later whether neutrinos are Majo- 
rana or Dirac particles 6, 7, 8, 9, 10]. Moreover, the rate 
of the Ov(3f3 decay, or Hmits on it, can tell us about the 
absolute neutrino-mass scale and to some extent about 
the neutrino mass hierarchy^ . But to achieve these goals 
we need an accurate evaluation of the nuclear matrix el- 
ements that govern the decay. 

In this paper, which builds on previous publications 
[T^ . [T^ (which we call I and II), we analyze some of 
the physics affecting the nuclear matrix element M^"^ 
— the competition between pairing and neutron-proton 
particle-particle correlations, the nonintuitive depen- 
dence of the decay amplitude on internucleon distance, 
and the treatment of short-range correlations and other 
high-momentum phenomena — that have not been suf- 
ficiently discussed before. As in our earlier papers (and 
most attempts to evaluate M'^'^) we use the Quasiparti- 
cle Random Phase Approximation (QRPA) and its gen- 
eralization, the Renormalized QRPA (RQRPA), with an 
interaction obtained from the G matrix associated with 
the realistic Bonn CD nucleon-nucleon interaction. That 
interaction, slightly renormalized, is used both as the like 
particle pairing and as the neutron-proton force. Where 
appropriate, we compare the results to those of the com- 
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plementary Large-Scale Shell Model (LSSM). 

The paper is organized as follows: In the next sec- 
tion, after briefly summarizing the relevant formalism, 
we show that the final value of Af reflects two compet- 
ing forces: the like particle pairing interaction that leads 
to the smearing of Fermi levels and the residual neutron- 
proton interaction that, through ground state correla- 
tions, admixes "broken-pair" (higher-seniority) states. A 
partial cancellation between these interactions increases 
the sensitivity to their strengths. The same tendencies 
are present in the LSSM, as a recent paper shows p^ . 
(During the processing of this manuscript, a new paper 
on the LSSM appeared [T^, emphasizing the competi- 
tion again.) In section III we discuss the dependence 
of Al'^'^ on the distance between the two neutrons that 
are converted into two protons. We show that the com- 
petition mentioned above implies that only internucleon 
distances r^j < 2-3 fm contribute. That fact, not rec- 
ognized before, explains the sensitivity of the decay rate 
to higher order terms in nucleon currents, nucleon form 
factors, and short-range nucleon-nucleon repulsion. We 
show that the surprising dependence on the internucleon 
distance occurs not only in the QRPA but also in an ex- 
actly solvable model [rZi] that contains many ingredients 
of real nuclear systems. Short-range correlations have 
recently inspired a lively discussion ^19i . .2(1 . ^] and we 
devote Section IV to various ways of treating them. In 
Section V we present numerical results for nuclei of ex- 
perimental interest, that include a comprehensive anal- 
ysis, within the QRPA method and its generalization, 
of the total uncertainty of the 0^/3/? nuclear matrix el- 
ements, and compare with results of the LSSM. Section 
VI summarizes our findings. Finally, in Appendix A we 
present formulae for two ways of evaluating the matrix 
elements, one via the evaluation of unsymmetrized two- 
body matrix elements (the procedure usually used) and 
another one through the product of two one-body matrix 
elements. And in Appendix B we show how to calculate 



shell-model particle-hole decompositions so they can be 
compared with those calculated in the QRPA. 



II. FORMALISM AND MULTIPOLE 
DECOMPOSITIONS 

Throughout we assume that the Oi^/3/3 decay, if ob- 
served, is caused by the exchange of the Majorana neutri- 
nos, the same particles observed to oscillate. The half-life 
of the decay is then 
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where G^'^ {Eq , Z) is a precisely calculable phase-space 
factor and M'~"^ is the nuclear matrix element. The ef- 
fective Majorana neutrino mass {mpji) is related to the 
absolute mass scale and oscillation parameters through 
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where Ud is the first row of the neutrino mixing ma- 
trix and the and ai are unknown Majorana phases. Any 
uncertainty in M'^'^ makes the value of {rapp) equally 
uncertain. 

As stated above, we use the QRPA and RQRPA meth- 
ods based on the G matrix derived from the realistic 
Bonn CD nucleon-nucleon force, i.e., the many body 
hamiltonian is 
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We describe in detail in Section V below the input used to 
solve the corresponding well known equations of motion. 

In the QRPA (and RQRPA) M^" is written as a sum 
over the virtual intermediate states, labeled by their an- 
gular momentum and parity and indices ki and kf 
(explanations of the notation are in Appendix A, and 
II): 
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FIG. 1: Contributions of different angular momenta J' asso- 
ciated with the two decaying neutrons to the Gamow- Teller 
part of M"" in *^Se (upper panel) and ^^"^Te (lower panel). 
The results of LSSM (dark histogram) 23] and QRPA treat- 
ments (lighter histogram) are compared. Both calculations 
use the same single-particle spaces: (/5/2,P3/2,Pi/2, <79/2) for 
^^Se and (37/2, ^5/2, ^3/2, Si/2, ^11/2) for "'^Te. In the QRPA 
calculation the particle-particle interaction was adjusted to 
reproduce the experimental 2;//3/3-decay rate. 
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The operators Ok,K = Fermi (F), Gamow-Teller (GT), 
and Tensor (T) contain neutrino potentials and spin and 

k ■ ' 

isospin operators, and RPA energies Ejl^ 
trino potentials, in turn, are integrals over the exchanged 
momentum q, 

HK{ri2,E'}^)^ 
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The functions fF,GT{qr 12) = joiqru) and fxiqr 12) = 
j2iqri2) are spherical Bessel functions (the sign of ^2 was 
given incorrectly in Ref. [l3|)- The functions hxiq^) are 
(4) defined in Appendix A and in II. The potentials depend 
explicitly, though rather weakly, on the energies of the 
virtual intermediate states, Ej^. The function ,f{ri2) in 
Eq. 0] represents the effects of short range correlations. 
These will be discussed in detail in Section IV. 

Two separate multipole decompositions are built into 
Eq. (|4]). One, already mentioned, is in terms the of 
the virtual states in the intermediate nucleus, the good 
quantum numbers of the QRPA and RQRPA. The other 
decomposition is based on the angular momenta and par- 
The neu- ities J^^ of the pairs of neutrons that are transformed into 
protons with the same J^'^ (we drop the superscript tt 
from now on for convenience). This latter representation 
is particularly revealing. In Fig. [T]we illustrate it both in 
the LSSM and QRPA, with the same single-single parti- 
cle spaces in each. These two rather different approaches 
{Ei + Ef)/2 ' agree in a semiquantitative way, but the LSSM entries 
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FIG. 2: Contributions of different angular momenta J asso- 
ciated with the two decaying neutrons to M^"^ in ^^Ge, ^""Mo 
and ^'^'^Te. We use the QRPA, with the interaction strength 
Ppp adjusted so that the 2vfij5 lifetime is correctly reproduced. 
Short-range correlations are included the same way as in I and 
II. 



for > are systematically smaller in absolute value. 

Ref. [IBl makes the claim that QRPA results are too 
large because they omit configurations with seniority 
greater than 4, which are especially effective in cancel- 
ing the pairing part of the matrix element. This state- 
ment is not correct. The QRPA does include configu- 
rations with higher seniority (4,8,12, etc.) and, as the 
Fig. [1] shows, the broken pair contributions to the ma- 
trix elements are as large or larger than in the LSSM. 
(Some of the difference might be due to differences in 
single-particle energies and occupation numbers, which 
are not identical in the two calculations even though the 
single-particle wave functions are.) The reason that the 
QRPA results presented in Section V are larger are a 
somewhat greater pairing contribution and contributions 
from negative parity multipoles that reinforce it. Most 
of the negative-parity contributions are absent from the 
shell model because of restrictions on the model space. 
These results suggest that the shell model is as likely to 
be missing important physics as is the QRPA. We return 
to this point in Section V. 

In Fig. Owe show the J decomposition for three nuclei 
in the QRPA, with single-particle spaces encompassing 
two major shells, a more natural span for this method. 
The cancellation between components with J = Q pairs 
and with J ^ ^ pairs is always pronounced. The net M'^^ 
is considerably smaller than the pairing contribution and 
so depends rather sensitively on the pairing physics that 
determines the J = Q part, as well as on the strength 
of the proton-neutron force that determines the J ^ 
part. 

From the structure of Eq. ^ and in particular from the 
form of the reduced matrix elements (Oj^|| [Cp,c„'] j| | J'^/c/) 
and (J'^fci 1 1 [CpC„] J 1 10^) it is obvious that only one 
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FIG. 3: The contributions of different intermediate-state an- 
gular momenta J to M"'' in ^""Mo (positive parities in the 
upper panel and negative parities in the lower one). We show 
the results for several values of gpp. The contribution of the 1"*" 
multipole changes rapidly with gpp, while those of the other 
multipoles change slowly. 



of the two possible couplings between the neutron 
and proton operators in the two-body matrix element 
{p{l),p'{2)- J II /(ri2)Oif/(ri2) II n(l),n'(2);^) is re- 
alized. This means that this two-body matrix element 
should not be antisymmetrized. In the LSSM one typ- 
ically uses the closure approximation, which represents 
Af^" as the ground-state-to-ground-state transition ma- 
trix element of a two-body operator. M^'^ can then be 
rewritten purely in terms of the antisymmetrized two- 
body matrix elements. After antisymmetrization, how- 
ever, it is not possible to recover the decomposition into 
the multipoles J'^ of the virtual intermediate states. In 
Appendix B we show how shell- model practitioners, by 
retaining unsymmetrized matrix elements, can decom- 
pose the matrix element into intermediate-state multi- 
poles J'^ for comparison with QRPA calculations. We 
cannot, however, make the comparison here without 
more shell-model data than has been published. 

When using the QRPA or RQRPA to evaluate M^" , 
one must fix several important parameters, the effects of 
which were discussed in detail in I and II. The strength 
Qpp by which we renormalize the Bonn-CD G matrix in 
the neutron-proton particle-particle channel is particu- 
larly important. We argued in I and II that gpp should 
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FIG. 4: The dependence on na of M"'' for ^"Ge, '""Mo and 
isorpg rpj^g upper panel shows the full matrix element, and 
the lower panel shows separately 'pairing' {J = for the two 
decaying neutrons) and 'broken pair' {J 7^ 0) contributions. 
The integrated matrix element is 5.35 for ^®Ge, 4.46 for '""Mo, 
and 4.09 for '^°Te. The Qpp values that reproduce the known 
Tj^/j are 1.030, 1.096 and 0.994. The single-particle space for 
^^Ge contains 9 levels (oscillator shells A'^ = 3, 4), and that for 
'°°Mo and '^°Te contains 13 levels (oscillator shells = 3, 4 
plus the / and h orbits from = 5). Short-range correlation 
are not included, i.e. f{ri2) = 1 in Eq. Q. 



be chosen to reproduce the rate of two-neutrino f3(3 decay. 
This choice, among other things, essentially removes the 
dependence of M'^'^ on the number of the single-particle 
states (or oscillator shells) in the calculations. The 
matrix element depends only on the 1"*" multipole. In 
Fig. [3] we show that it is essentially this multipole that 
is responsible for the rapid variation of Al'^'^ with gpp. 
Fixing its contribution to a related observable {2i^ de- 
cay) involving the same initial and final nuclear states 
appears to be an optimal procedure for determining gpp. 



III. DEPENDENCE ON THE DISTANCE 
BETWEEN THE NUCLEONS INVOLVED IN 
THE OuPP TRANSITION. 

The operators Ok in Eq. dH) depend on the distance 
ri2 between the two neutrons that are transformed into 
protons. The corresponding neutrino potentials arc the 
Fourier transforms over the neutrino momentum q as 




FIG. 5: The ri2 dependence of M^" for ^"^Ge, from calcu- 
lations with different number of single-particle orbits. The 
dot-dashed curve was obtained with 21 s.p. subshells, the full 
curve with 12 subshells, the dashed curve with 9 subshells, 
the dotted curve with 6 subshells and the double dot-dashed 
curve with only 4 subshells. 



shown in Eq.(I5]). Obviously, the range of ri2 is restricted 
from above by ri2 < 2Rnuci- We show here, however, 
that in reality only much smaller values, ri2 ^ 2-3 fm, or 
equivalently larger values of q, are relevant. Thus a good 
description of the physics involving distances ri2 ~ 1 fm, 
or q ^ 200 MeV is important. That finding has not been 
recognized before, but perhaps it should be not so sur- 
prising that q PFermi is the most relevant momentum 
transfer. 

An example of the ri2 dependence of M'^'^ is shown in 
Fig. [4] for three nuclei. The quantity C(r) is defined by 
evaluating M"" after multiplying H^{t', Ej^ ) by r'^S{r — 
r'), so that C(r) is the contribution at r to Al'^" , with 
C{r)dr = AI^'^ . As the lower panel of the figure 
demonstrates, the cancellation between the — and 
J ^ ^ components is essentially complete for r\2 ^ 2-3 
fm. Since the typical distance from a particular nucleoli 
to its nearest neighbor is ~ 1.7 fm (because Rnuci — 
1.2A^/'^) the nucleons participating in the Oi'PP decay are 
mostly nearest neighbors. Short-range nucleon-nucleon 
repulsion, the finite nucleon size, represented by nucleon 
form factors, and components of the weak currents that 
are typically suppressed by q/Mnucieon are therefore more 
important than one would naively expect. 

Perhaps the most interesting thing about the figure is 
that the pairing and non-pairing parts of C(r) taken in- 
dividually (as in the two panels of the figure) extend to 
significantly larger r. The cancellation between them, 
that we discussed earlier, is particularly effective beyond 
2 or 3 fm, leaving essentially nothing there. Figure O 
shows that the shape of C(r), like the integrated ma- 
trix element, is essentially independent of the number of 
single-particle orbits included, as long as the truncation 
is not too severe (as it is with the dash-double-dot curve, 
for which important spin-orbit partners were omitted - 
only the 4 single particle states ^3/2,^1/2, /5/2, 59/2 were 
included) and the coupling constant gpp is chosen to re- 



FIG. 6: The ri2 dependence of the contributions of various 
pieces of M°'' for ^^Ge are shown. Here AA stands for axial, 
VV for vector, AP for axial-pseudoscalar interference, PP for 
pseudoscalar, and MM for weak-magnetism. For definitions 
see Appendix A, or Ref. [2^. The model space contains 9 
subshells. 

produce the measured 2i^/3(3 lifetime. For other values of 
gpp the cancellation between the — and J ^ ^ con- 
tributions at r larger than 2 or 3 fm is not as complete 
as in Fig. U) We return to this point shortly 

We show the r\i dependence of the different parts of 
the M^^ in Fig. [6l All individual contributions die out 
at r larger than 2 or 3 fm. The pseudoscalar-axial vector 
interference part has opposite sign from the other con- 
tributions, and essentially (and accidentally) cancels the 
contributions of the vector, weak magnetism and pure 
pseudoscalar pieces. The higher-order terms reduce the 
matrix element noticeably, and have to be included. 

To gain some insight into the renormalization of the 
double-beta decay operator in the shell model, Ref. [T7| 
employs a solvable model based on the algebra 50(5) x 
50(5). The valence space contains two major shells 
ifP99/2 a-nd 5^37/2), split by an energy e, with degenerate 
levels within each shell. The single-particle wave func- 
tions are taken from a harmonic oscillator with huj =9.2 
MeV. The Hamiltonian (for this schematic model only, 
not in the rest of the paper) is 

H = €N2 (6) 
2 

~G ^ ^ i^pp^pp ~^ ^nn^nn ~^ dpp^pn^pn ~ 9phTa ' Tf,) , 

a,b=l 

where a,b = 1,2 label the shells (lower and upper), e is 
the energy difference between the shells, N2 is the number 
operator for the upper shell, Ta is total isospin operator 
for shell a, and 

PP ~ 2 / ^ U'a^alO 
afEa 
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FIG. 7: The ri2 dependence of M"" in the exactly solvable 
model for 4 values of gpp (upper panel). The integrated matrix 
elements are 2.93 for gpp = and 1.69 for gpp = 1. The lower 
panel shows separately the contributions of J — pairs and 
J" / pairs for gpp = 0.8. 

4n = ^E^-t^Hir (7) 

Here tt^ (i^l,) creates a proton (neutron) in level a with 
angular momentum ja, j = ^/2j + 1, and the square 
brackets indicate angular-momentum coupling. H con- 
tains only generators of SO{5) x 50(5), so its lowest lying 
eigenstates consist of configurations in which the nucle- 
ons are entirely bound in isovector 5 pairs of the type in 
Eq. 0. 

This model has no active spin, so it is only suitable 
for calculating Fermi (neutrinoless) double-beta decay. 
To simulate the effect of gpp on Gamow- Teller decay we 
change the Fermi matrix element by varying the strength 
of isovector neutron-proton pairing, in the same way that 
we change the Gamow- Teller matrix element in the real- 
istic QRPA by varying the strength of isoscalar pairing. 
The advantage of this model is that we can solve it ex- 
actly rather than in the QRPA. 

We can use the model to test the ri2-dependence of the 
double-beta decay matrix element in an exact solution. 
(Analytic expressions for the necessary matrix elements 
are in Ref. [i3|-) The upper panel of Fig. [7] shows the 
dependence for several values of gpp, with gph = and 
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FIG. 8: The momentum-transfer dependence of M in Ge. 
The upper panel is for the full matrix element; in the lower 
panel we separate the J = and J parts. The scale 
is different in the two panels. The model space contains 9 
subshells. 



e = lOG. Just as in the realistic QRPA calculations, the 
contribution beyond r = 3 fm is very small for gpp around 
1; it is too small to distinguish from zero in the figure be- 
yond 5 fm for gpp = 0.8. However, for other values, as 
noted above, the large-r contributions can be substantial. 
The bottom panel divides the function into like-particle 
pairing and non-pairing parts for gpp = 0.8. The two 
cancel to high precision at large r. The suppression at 
long ranges we observe in the QRPA, then, appears to be 
fairly general. It happens even in a very simple model, 
solved exactly. (As noted above, a new preprint [l^ ap- 
peared during the processing of this manuscript. In it the 
ri2 dependence of the M'^'^, as well as the dependence of 
the separated pairing and broken-pair contributions, was 
evaluated in the LSSM. That analysis, inspired by our 
work, yielded curves that are strikingly similar to those 
in Fig. 12) 
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FIG. 9: The ri2 dependence of M°'" in '^'^Ge evaluated in the 
model space that contains 9 subshells. The four curves show 
the effects of different treatments of short-range correlations. 
The resulting M"" values are 5.32 when the effect of short 
range correlations is ignored, 5.01 when the UCOM trans- 
formation [2J is applied, 4.14 when the J{r 12) from Fermi 
hypernetted-chain calculations [2^ is used in Eq.Q, and 3.98 
when the phenomenological Jastrow f{r 12) is used [26| . 



IV. SHORT-RANGE CORRELATIONS AND 
OTHER HIGH-MOMENTUM PHENOMENA 

Since only ri2 <i 2 — 3 fm, (equivalently q > hc/{2 — 
3 fm)), contributes to M""^, some otherwise negligible ef- 
fects become important. These effects are not commonly 
included, or included only in rough approximation, in 
nuclear-structure calculations. For example, the dipole 
approximation for nucleon form factors and the corre- 
sponding parameters My and Ma come from electron- 
and neutrino charged-current-scattering from on-shell 
nucleons. Nuclear structure deals with bound nucleons 
and virtual neutrinos that are far off-shell. Similarly, the 
induced pseudoscalar current, with its strength obtained 
from the Goldberger-Treiman relation, has been tested 
in muon capture on simple systems. Here we are using 
this current for off-shell virtual neutrinos. Short range 
nucleon-nuclcon repulsion has been considered carefully 
when calculating nuclear binding energies, but here we 
need its effect on a transition operator connecting two 
different nuclear ground states. All these effects will in- 
troduce some uncertainty because their treatment is not 
well tested. Nevertheless, it is important to understand 
their size at least roughly. 

To show the importance of high momenta explicitly, 
we display in Fig. [5] the q dependence C{q) — defined in 
complete analogy to C(r) — of M°'^ in ^^Ge, in a similar 
manner to which we exhibited the ri2 dependence earlier. 
The cancellation between the = and J parts is 
particularly complete at lower values of q so that the 
resulting curve in the upper panel, although reduced in 
magnitude, is clearly shifted towards higher q. 

The first high-momentum effect we examine is short- 
range correlations. Fig. [5] displays the rn dependence of 
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FIG. 10: The dependence of IVr" in ^"Ge on the value of 
the dipole-form-factor cut-off parameters Ma,v- The upper 
curve was calculated without short range correlations, while 
the two lower curves were obtained with the UCOM and Jas- 
trow methods. The two lower curves, unlike the upper one, 
are essentially flat for Ma,v larger than the standard value of 
^ 1 GeV. 
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FIG. 11: Circles represent the mean value of the upper and 
lower limits of our calculated values of M see text for 
description of error bars. For comparison the results of a 
recent Large Scale Shell Model evaluation of M that used 
the Jastrow-type treatment of short range correlations are 
also shown as triangles 



Af '^^ for several methods of handling short-range physics. 
For obvious reasons aU methods reduce the magnitude 
of M^". The Unitary Correlation Operator Method 
(UCOM) leads to the smallest reduction, less than 
5%. The phenomenological Jastrow-like function /{ru) 
in Eq.dll) (from Ref. (H) reduces M^" by about 20%. 
We also display the results of using a microscopically- 
derived Jastrow function [2^ ; its effect is similar to that 
of the phenomenological function. Since it is not clear 
which approach is best, we believe it prudent to treat 
the differences as a relatively modest uncertainty. 

Nucleon form factors pose fewer problems because it 
turns out that once the short-range correlations effects 
are included, no matter how, the form factors are almost 
irrelevant as long as the cut-off masses Ma.v are at least 
as large as the standard values (Ma = 1-09 GeV and 
My = 0.85 GeV). In Fig. [10] we show the dependence 
of M^^ on the values of Ma,v which for this purpose 
are set equal to each other, with three alternatives for 
treating short-range correlations. By 2 GeV the curves 
have essentially reached the infinite-mass limit. Since 
they are essentially flat past 1 GeV for both the UCOM 
and Jastrow-like prescriptions, including the form factors 
causes only minor changes in Af"''. Only if the correla- 
tions are ignored altogether do the form factors make a 
significant difference. 

Finally, there is little doubt that the higher order weak 
currents, induced pseudoscalar and weak magnetism, 
should be included in the calculation. Even though the 
Goldberger-Treiman relation has not been tested in two- 
body operators, the relation is sufficiently well estab- 
lished that we do not associate a sizable uncertainty with 
its use. 



V. NUMERICAL RESULTS 

Fig-HDshows our calculated ranges for Af'"'^, defined as 
AI'""" = (.gA/1.25)2Af°^ to allow us to display the effects 
of uncertainties in gA- Such a definition allows us to use 
the same phase space factor G^'^ with qa = 1-25 when 
calculating the Oi'/S/J-decay rate. 

We consider the effects of short-range correlations an 
uncertainty, since we don't know the best way of treat- 
ing them. The error bars in Fig. [11] represent the dif- 
ference between the highest and lowest of 24 calculations 
— in either the QRPA or RQRPA, each with 3 different 
sets of single-particle states (usually 2, 3 and 4 oscillator 
shells), 2 values for the in-medium gA (1-0 and 1.25), and 
2 treatments of short-range correlations (phenomeologi- 
cal Jastrow functions and the UCOM method) — and 
include the experimental uncertainty in the values of the 
2i> lifetimes used to determine gpp. Thus the error bars 
displayed in Fig. [TT]represent our estimate of the full un- 
certainty in the Qvj3[3 matrix elements within the QRPA 
and RQRPA methods. Though the results are in reason- 
able agreement with those presented in Refs. [lO, [2l|, the 
uncertainties are different. 



Reference 27[ presents new LSSM results for Af'"'^, 
the values of which are shown also in Fig. [11] They are 
somewhat smaller than the QRPA values for "^^Ge and 
^^Se and in a fair agreement for "''Cd, ^^Sxe, ^'^'^Te and 
^^^Xe. Various nuclear-structure effects may be respon- 
sible for the discrepancies. More complicated configura- 
tions in LSSM that are absent in QRPA typically reduce 
Ai^^ and could make most of the difference. On the other 
hand, in the QRPA includes more single particle states 
than in LSSM. That has the tendency to increase M^^ 
and could also be responsible for the discrepancy. As 
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we said earlier, the notion that the QRPA omits high- 
seniority states is not correct and shouldn't be used to 
argue that the LSSM results are more accurate. 

Given the interest in the subject, we show the range 
of predicted half-lives corresponding to our full range 
of M'O" in Table U (for {mpp) = 50 meV). As we ar- 
gued above, this is a rather conservative range within 
the QRPA and its related frameworks. One should keep 
in mind, however, the discrepancy between the QRPA 
and LSSM results as well as systematic effects that might 
elude either or both calculations. 



hospitality during the initial stage of this work and the 
US Department of Energy for partial support. 



APPENDIX A 

Here we outline the derivation of Af , emphasizing in- 
duced (higher-order) currents. We assume light-neutrino 
exchange throughout, and the standard lepton x hadron 
weak charged-current Hamiltonian. The hadronic cur- 
rent, expressed in terms of nucleon fields is 



VI. CONCLUSIONS 

The most important and novel result here is that the 
generic competition between — Q (pairing) and J ^ Q 
(broken pair) multipoles leads to almost complete can- 
cellation of the contribution to the matrix element from 
internucleon distances r > 2-3 fm. That explains why 
the effects that depend on smaller values of r, or equiva- 
lently larger momentum transfers g, become important. 
This competition also means that the final matrix ele- 
ments have enhanced sensitivity to the strengths of these 
interactions. Despite the uncertainties associated with 
the short range effects, we conclude that a proper fitting 
of the QRPA and/or RQRPA parameters leads to a rela- 
tively narrow range for M'^^ ^ with a smooth dependence 
(^^Zr being an exception) on atomic charge and mass. 

We evaluate the values of the matrix elements for nu- 
clei of experimental interest and display our best estimate 
of the corresponding spread. Part of that spread is asso- 
ciated with the difference in the size of the single particle 
space and whether QRPA or RQRPA is used, as dis- 
cussed earlier in I and II. An interesting new conclusion 
is that short-range correlations, no matter how they are 
treated, essentially eliminate the effect of finite nucleon 
size on the matrix elements. But we still do not know the 
best way to treat the correlations, a fact that contributes 
about 20% to uncertainties presented above. The uncer- 
tainty in the effective value of contributes about 30%, 
with the rest due to choice of method and model space, 
and the experimental uncertainty in 2i/ lifetimes. 



Acknowledgments 



5y(g )Y + i9M{q 
~9A{q^)Yl^ - 9p[q^Wl^] (Al) 



where rup is the nucleon mass and is the momen- 
tum transfer, i.e. the momentum of the virtual neu- 
trino. Since in the Qvfip decay (7^ ^ (/g we take 
g2 ^ —(f. For the vector and axial vector form fac- 
tors we adopt the usual dipole approximation gviq"^) = 
gv/{l + q^/M^)\ gA{q^) = W(l + 9 7^1)2, with 
gy = 1, = 1.254, Mv = 850 MeV, and Ma = 
1086 MeV. We use the usual form for weak magnetism, 
gM[q'^) = {fip—fin)gv{q^), and the Goldberger-Treiman 
relation, gp{q^) — 2'mpg A{q / {q + m^) , for the in- 
duced pseudoscalar term. 

To derive the expression for the matrix element we 
follow the procedure outlined in Ref. [1^, arriving after 
a few steps at an expression for the 0^ 0^ ground 
state to ground state transition: 



At 



AttR 
IT 



1 



(27r)3 



-iq.(xi-X2) 



(A2) 



< 0+1 Jt (fi)|m X m| J"t(f2)|0+ > 
? E^-iE.^E,)/2 + \q\ ^'^^^^^^^ ■ 

We have made the (accurate) approximation that all elec- 
trons are emitted in the Si/2 state, with energies equal 
to half the available energy {Ei — Ef)/2. The normal- 
ization factor 47r R/g\, introduced for convenience, is 
compensated for by corresponding factors in the phase 
space integral. 

Reducing the nucleon current to the non-relativistic 
form yields (see Ref.^]) in Eq.(|M|: 
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(A3) 



n— 1 k 



where J"{q ) — gv{q ) and 



jf^2\ i^2\-'^n^q I f^2\^ /->2\q^ri'q 

Jn[q ) = 9M{q n-z V9A{q )cr ~ gp{q )—z , 

Znip Znip 

(A4) 

r,i is the coordinate of the nth nucleon, k — 1,2,3, and 
gP " is the metric tensor. 
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TABLE L The calculated ranges of the nuclear matrix element M '''' evaluated within both the QRPA and RQRPA and with 
both standard {qa ~ 1.254) and quenched (gA ~ 1-0) axial-vector couplings. In each case we adjusted gpp so that the rate of 
the 2i^/3/3=decay is reproduced. Column 2 contains the ranges of M "'^ with the phenomenological Jastrow-type treatment of 
short range correlations (see I and II), while column 3 shows the UCOM-based results (see Ref. [13]). Columns 3 and 5 give 
the 0!//3/3-decay half-life ranges corresponding to the matrix-element ranges in columns 2 and 4, for < mpp >= 50 meV. 

Nuclear (R)QRPA (Jastrow s.r.c.) (R)QRPA (UCOM s.r.c.) 



transition M {{mpp) ^ 50 meV) M {{mpp) ^ 50 meV) 

'"Ge '"Se (3.33,4.68) (6.01, 11.9) x 10'"' (3.92,5.73) (4.01,8.57) x 10^" 

^'^Se^^^Kr (2.82,4.17) (1.71, 3.73) x lO^*^ (3.35,5.09) (1.14, 2.64) x 10^'' 

96 2r ^^Mo (1.01, 1.34) (7.90, 13.9) x 10^'^ (1.31, 1.79) (4.43, 8.27) x 10^'' 

^°°Mo^ i°°i?u (2.22,3.53) (1.46,3.70) x 10^*^ (2.77,4.58) (8.69,23.8) x 10^^ 

"'^Cd^^^S'n (1.83,2.93) (1.95, 5.01) x lO^'' (2.18,3.54) (1.34, 3.53) x 10^*^ 

i28yg^i28j^g (2.46,3.77) (3.33, 7.81) X lO^'^ (3.06,4.76) (2.09, 5.05) x lO^'' 

i30yg^i30j^g (2.27,3.38) (1.65, 3.66) x 10^*^ (2.84,4.26) (1.04, 2.34) x 10^*^ 

"^Xe ^ ^^^^Ba (1.17, 2.22) (3.59, 12.9) x 10^^ (1.49,2.76) (2.32, 7.96) x 10^^ 



The two current operators in M'^'^ lead to an expression 
in terms of 5 parts \22\ : 



M""" = Mvv + M 



Maa + Map + Mpp, (A5) 



with the notation indicating which parts (axial, vector, 
etc.) of the nucleon current contribute. After integrating 
over dxi, dx2 and dftg in (jA3[) and writing one-body 
charge-changing operators in second quantization as 



,A V- (pII H ^) r 1 



(A6) 



we obtain 



J,7r,ki ,kf pnp' n' 



(A7) 



R 

9a Jo 



, , j(q) 

pnp n , J V^/ 



hK{q^)q^dq x 



mm -r 

{o^\\[4£n']j\\J^kf){rkf\rk,)(rk,\\[4~c,,]j\ 

Here K = VV, MM, AA, PP, AP and 



hvv{q^) = -gviq"^)-, 



I'MM 



hAA{q') = 9\{q^), hpp{q^) = 

,9A{q^)gp{q'^)q^ 



9M{q'^)q'^ 

gp{q^)q'^ 
^rrip 



hAp{qn = -2- 



2to„ 



(A8) 



The reduced matrix elements of the one-body opera- 
tors c+c„ (the tilde denotes a time-reversed state) in Eq. 
(|A8|) depend on the BCS coefficients Ui,Vj and on the 
QRPA vectors X, Y [2j| . The nuclear structure informa- 
tion resides in these quantities. 

The 'Pp^p,,^, j{q) in Eq. ((X8|l are products of the re- 
duced one-body matrix elements of operators ©'^"^(g): 

-P^Z'n'Al) = (P II II ^)(P' II Oj^il) II 



^ pnp' n' ^ J ^ ^ ( ' 



J+L + 1 



X 



L=J,J±1 

{p\\ Of](g) \\n){p' II 0[%) II n'), 
-PZ'n'jiQ) = {P II Of{q) II n){p' II Of\q) \\ n'), 
-p^P, , Jn) = r^^, , Aq) 

' pnp'n'..J\^J ' pnp' n' .J f 

'pnp'n'.J\HJ ' pnp' n' .Jy^J ' pnp'n' ,J\^J ' 



(A9) 



Here 



0%{q) = 2V2 j J {qr)YjM{nr), 
Oflj,j{q) = 2V2jLiqr){YL{nr)®a,}jM, 

0f],,(g) = 2V2^l^^jj^,{qr)Cj\, x 

{Yj^liflr) «) (Ji}jM 
I2 J -\- 

-2V2^^j^JJ+i{qr)CZio 10 X 

{rj+i(f7,)®ai}.7A/- (AlO) 

The final step, leading to Eq. ((4]) in the text, is to 
rewrite the product of two one-body matrix elements 
as an appropriately recoupled (with pairs of protons 
and neutrons coupled to angular momentum J^) unsym- 
metrized two-body matrix element. Without the compli- 
cations of angular momentum, this step simply reads 

(p|0(l)|n)(p'|0(2)|n') = (p,p'|0'(l, 2)|n, n'). (All) 

We then transform to relative and center-of-mass coor- 
dinates ri2 = ri — r2 and R12 = {ri + r2)/2. Since the 
exchange potential depends only on ri2 = |fi2| we end 
up with Eq. g]). The Fermi (F), Gamow-Teller (GT) and 
Tensor (T) operators in that equation are 

0F{n2,E%) = T+{l)T+{2)HF{n2,E'j^), (A12) 
OGT{ri2,E%) = T+{l)T+{2)HGT{ri2,E'}^)ai2, 
OT(ri2,^^) = T+(l)T+(2)i7T(ri2,^^)5i2. 
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Here 

0'12 = f?l • CT2 , 

5l2 = 3((Ti • fi2)(l72 • ri2) - (Ti2 . (Al3) 

The functions hK{q^) that determine the Hk's through 
the integrals over q in Eq. ([5]) are 

hpiq'^) = -9Ahvv{q^), 
hcriq^) = ^{2hMM{q^) + hpp{q^) + hAp{q^)) 
+hAAiq^), 

hriq^) = -^{hMMiq^) - hpp{q^) - hApiq^)) , 



and the fuU matrix element is 



Mp 



M"" = 4-+ Mgt + Mt 

9 a 



(A14) 



(A15) 



Short range repulsion can then be included as explained 
in Section IV. 



APPENDIX B 

Here we show how to calculate shell-model particle- 
hole decompositions so they can be compared with those 
calculated in the QRPA. To avoid too many complica- 
tions, we will use the closure approximation. (In the 
text we have shown that, within the QRPA at least, us- 
ing the closure approximation for the Oz^/3/3-decay results 
in an error of < 10%.) The matrix element M (the sub- 
script K is implied) can be written as in Eq. gl), with the 
overlap between intermediate-nucleus eigenstates a Kro- 
necker delta if the those states are determined uniquely 
(as in the shell model). The matrix element M can be 
decomposed: 



M 



(Bl) 



where a stands for the set of indices p,p' , n, n' and 

A^" = E4JO^ (B2) 



jj 



The parity index tt is implicitly included along with J 
and J . The O" are given by 



Y.{o^[4rcn-]j\\rkf){rkf\rh) 



ki,kf 



x(J^fc.||[c+£„]j||0+) , (B3) 
and the Sjj are everything else in Eq. (jH): 



with J EE ^2J + 1 and 

Z<} = {p{l),p'{2)-J II />i2)Oif/(ri2) II n(l),n'(2); J) . 

(B5) 

(The f{ri2) can be omitted if short-range correlations are 
included some other way.) 

Now we can write M" in two different ways: 



with 



and 



(B6) 



(B7) 



(B8) 



The Mj''" are the pp-hh amplitudes and the Mj'^'" are 
the ph multipole-multipole amplitudes that we want to 
calculate in the shell model. All the nuclear structure 
information is in the Oj . 



From Eq. (|B7|) we have 



J 



and 



(B9) 



(BIO) 



J 



(from Eq. (|B9 



So, exchanging the primed and unprimed labels in the 
sum, we can write the AI^j^'" in terms of the Mj''" as 



M- 



J,J' 



The final particlc-holc multipole contribution that we 
want is then just 



From the relation 



X ^ 



( a b 




fab 




\c d 




{c d 





(B12) 



Sj,j' (B13) 



we have 



] j\Z}, (B4) sjt^ 



, (B14) 
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Finally, putting everything together, we get 



M 



ph 



p,p' ,n,n' ,J J' 



(B15) 



In a shell model calculation, one can write the double 
beta-decay matrix element solely in terms of antisym- 
metrized matrix elements of the corresponding operator. 
But to obtain the particle-hole decomposition above, the 



natural definition since the operator really represents a 
product of two one-body currents, one must start from a 
representation of the operator in terms of unsymmetrized 



matrix elements Zj 



O 



K 



E 

p.n,p' ,n' 



ryOi 



1 



(B16) 



and calculate the M^'"" , for all a = p,p',n,n', not just 
p > p' , n > n' . 
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